arXiv:astro-ph/0211492vl 22 Nov 2002 


Submitted to ApJ. 


On the Eccentricity Distribution of Coalescing Black Hole 
Binaries Driven by the Kozai Mechanism in Globular Clusters 

Linqing Wen^ 

Division of Physics, Mathematics, and Astronomy 
California Institute of Technology 

February 2, 2008 

ABSTRACT 

In a globular cluster, hierarchical triple black hole systems can be produced 
through binary-binary interaction. It has been proposed recently that the Kozai 
mechanism could drive the inner binary of the triple system to merge before it 
is interrupted by interactions with other field stars. We investigate qualitatively 
and numerically the evolution of the eccentricities in these binaries under 
gravitational radiation (GR) reaction. We predict that ~ 30% of the systems 
will possess eccentricities >0.1 when their emitted gravitational waves pass 
through 10 Hz frequency. The implications for gravitational wave detection, 
especially the relevance to data analyses for broad-band laser-interferometer 
gravitational wave detectors, are discussed. 


Subject headings: gravitational waves—binaries: close—stellar dynamics— 
relativity 


1. Introduction 

Globular clusters are excellent birthplaces for black hole (BH) binaries. In our Galaxy, 
a globular cluster (GG) contains typically N = 10® stars and resides normally in the galactic 
halo. Typical ages of these clusters are estimated to be 12 billion years, and the clusters 
contain the oldest stars known in our Galaxy. Stellar evolution theory predicts that massive 
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stars with mass above 10 Mq would evolve into black holes or neutron stars in ~ 10^ years. 
For a typical GC that follows the initial mass function (IMF) given by Scalo (1986), a 
signihcant number of BHs (~ 6 x 10“^iV for initial masses > 20-25Mq) should have been 
born. These BHs would be the most massive stars in the cluster, and hence dynamical 
two-body relaxation causes them to sink rapidly towards the core. In the core, BH binaries 
are expected to form. The single BHs form binaries preferentially with another BH, and 
BHs born in binaries with a lower mass stellar component rapidly exchange the companion 
for another BH through three body relaxation processes. If these BH binaries merge within 
the Hubble time, they would become excellent sources for upcoming gravitational wave 
(GW) detectors (e.g. [Thorne 1987|) . 


It has been argued that globular clusters are unlikely places to harbor BH mergers 
(or subsequently to harbor a massive central black hole) if we consider binary-single 
interactions alone (e.g. Portegies Zwart fc McMillan 2000|) . Numerical simulations show 
that binary-single interactions tend to throw out these BH binaries and only 8% of BHs 
will be retained in the GG within its lifetime. Miller & Hamilton (2002) argued that this 
picture will be different if we consider binary-binary interactions. A substantial fraction of 
hierarchical triple BH systems can result from binary-binary interactions. The well-known 
Kozai mechanism ( jKozai 1962|) can then drive the inner binaries of a signihcant fraction 
of these triples to extreme eccentricities and thus reduce their time to coalescence to be 
shorter than the Hubble time. This is because the coalescence time scale tgr due to 
gravitational radiation for a binary with hxed semi-major axis a has a strong dependence 
on the eccentricity e through the quantity e = 1 — ( [Peters 1964]) , 


Tgr ^ 6 X 10 


10 


(a/AU)^(6/0.01)^-^ 

(mo -|- mi)mQmi/M q^ 


yr, 


( 1 ) 


where mo and mi are the masses of each binary component. In their paper. Miller & 
Hamilton (2002) have compared the coalescence time scale with the time scale for the 
system to be disrupted by interactions with held stars. They concluded that in a substantial 
fraction of cases, the system will achieve extreme eccentricities (e ~ 0) and will merge by 
gravitational radiation before the disruption. The BH mergers may thus be retained in the 
GGs and a central intermediate-mass BH may form through subsequent mergers. 


We are motivated by the fact that the proposed mechanism promises stellar mass 
BH-BH mergers with extremely high initial eccentricities and that they are potential 
sources for upcoming gravitational wave detectors such as LIGO. We investigate in this 
paper whether or not these systems would still possess hnite eccentricities when the emitted 
gravitational waves enter LIGO’s frequency band. The eccentricities of these sources are 
particularly interesting as no astrophysical merger sources were expected before to possess 
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significant eccentricities when they reach the LIGO band. Cnrrent efforts in detecting GWs 
from coalescing BH binary sources thus focus on waves from circular orbits (see |Buonanno( 
Ghen, fc Vallisneri 2UU2| for recent development). On the other hand, it is well-known that 
the gravitational waveform of an eccentric binary could be signihcantly different from that 
of a circular one and different detection strategies are required for optimal detections (e.g., 


Peters fc Mathews 1963| ; [Peters 1964|; [Lincoln fc Will 1990| ; [Moreno-Garrido, Mediavilla, fc 
Buitrago 1995[ ; [Gopakumar &: Iyer 2002[; [Glampedakis, Hughes, fc Kennehck 2002[) . Our 
work can thus help shed light on the importance of the eccentricities for detecting GWs 
from these sources. 


In this paper, we investigate in detail the eccentricity evolution of these merger 
systems under the impact of gravitational radiation reaction (simplihed as the GR effect). 
This is in addition to the Kozai mechanism and the post-Newtonian periastron precession 
(simplihed as the PN effect) considered by previous authors. It has been long realized 
( Peters &: Mathews 1963[ ; Peters 1964[) that gravitational radiation reaction is very efficient 
in circularizing coalescing binaries. For the systems we are considering, this circularization 
effect would have to compete with the “eccentricity enhancing” Kozai mechanism. Under 
the Kozai mechanism and the PN effect, the system tries to maintain secular cyclic 
eccentricity oscillations at hxed amplitude (with exceptions at critical inclination angles, 
see text). The GR effect operates most efficiently near the minimum e of each cycle. In 
some cases, it could drive a system near its minimum e directly towards merger. In other 
cases, it would affect the minimum e and alter the otherwise fairly strict periodic evolution 
trajectory of the system. In particular, it will enhance rapidly the importance of the PN 
effect which in turn would introduce fast oscillations to destroy the Kozai resonance. The 
Kozai resonance will eventually be quenched either directly by the GR effect, or by the PN 
effect before the GR effect takes over. We will discuss the details of these effects and their 
implications for the evolution of the eccentricities and gravitational wave frequencies. 


In section we give a summary of the relevant features of the Kozai mechanism and 
derive an explicit general expression for the minimum e (or maximum eccentricity) the 
system can reach for arbitrary values of BH masses and initial eccentricities. In section ^ 
we discuss the evolution of the system parameters under gravitational radiation reaction 
(the GR effect). The equations of motion are described in subsection We describe 
technical details to estimate the minimum e the system can reach in subsection [3.2[ . The 
PN effect and its role in determining the evolution of the gravitational wave frequency 
are described in subsection |3T3[ . A procedure to obtain the relation of the gravitational 
wave frequencies and the eccentricity at various evolution stages is given in subsection [3.4 
Two numerical examples are presented in subsection [3.5[ - In section we estimate the 
eccentricity distribution when the gravitational waves enters the LIGO band. Our results 
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are 


then summarized in section 


2. Kozai Mechanism 


2.1. Principle and Notation 


The Kozai Mechanism was initially proposed (P:4ozai 1962|) to describe the perturbation 
of asteroids due to Jupiter. This mechanism describes the secular evolution of the orbit of 
a hierarchical triple system as a result of cyclic exchange of angular momentum between 
the inner binary and orbit of the distant tertiary. One of the most interesting results of the 
Kozai mechanism is that the eccentricity of the inner binary can reach very large values 
while the semi-major axis remains unchanged, provided that the initial mutual inclination 
angle of the two binaries reaches a certain critical value. Although the mechanism was 
proposed for a system where the mass of one component in the inner binary is much smaller 
than those of the others, the mechanism has then been found to operate in hierarchical 
triple systems with arbitrary masses (|Lidov fc Ziglin 19761). 


A hierarchical triple star system consists of a close binary of masses mo, mi, and a 
tertiary m 2 in a much wider outer orbit (Fig. [^). The dynamics of the system can be 
treated as if it consists of an inner binary of point masses mo and mi, and an outer binary 
of point masses m 2 and Mi = tuq + mi. The separations of the masses each binary are 
denoted ri and r 2 . The mutual inclination angle of the two binaries is denoted I. For the 
inner and outer binaries respectively, we denote the eccentricities ei and 62 , reduced masses 
Pi = momi/Mi,p 2 = m 2 Mi/(Mi -|- m 2 ), semimajor axes ai and 02 . We dehne gi as the 
argument of pericenter of the inner binary relative to the line of the descending node. 

Throughout this paper, we frequently use quantities e = 1 — and 62 = 02 y^l — e| to 
simplify our discussion. Angular momentum is normalized to have dimension i/Length by 
dividing by the quantity Li = piv^GMJ. We use notation oq, / 3 o, and 7 for the normalized 
magnitude of total angular momentum of the triple system, of the outer binary, and of the 
inner binary respectively. Our ao and /do are related to dimensionless quantities a and (3 
used in previous literature by a = ao/y^, and [3 = /do/ydq. It follows that 7 = ^Oic, 

(3q = P 2 /hiY^M 2 a 2 /Mi(l — e|), and that cto is related to (3q and oie through the following 
trigonometric relation 

( 2 ) 

2 /do-y/aie 

As we shall see, during the binary’s Kozai mechanism-induced evolution, oq, (3q, and Oi are 
conserved while e changes. It follows that ao = 13q is the critical condition for the system to 
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be able to evolve to e = 0 (ei = 1). This implies a critical initial I of 


Io = Ic 



( 3 ) 


where Cq is the initial e. For this solution to exist, /3 > ^/eQ/2 is required. 

The total Hamiltonian of the system can be written as if = Hi + H 2 + H' , where 
ffi and H 2 are the Hamiltonians for the inner and outer binaries respectively as if they 
were isolated binaries of point masses. H' is the perturbative Hamiltonian which includes 
perturbation of the inner binary by the presence of the tertiary m 2 and perturbation 
of the outer binary by the hnite size of the binary component Mi. For ri/r 2 -C 1 and 
(m 2 /Mi)(ri/r 2 )^ 1, the Hamiltonian H' can be expanded in powers of ri/r 2 - The hrst 

non-zero contribution comes from the term of order (ri/r 2 )^, labeled as quadrupole term 
because of the degree of the Legendre polynomial associated with it. This term is sufficient 
for our discussion as we are only concerned about cases with large initial mutual inclination 
angle Jq (see |Ford, Kozinsky, fc Rasio 2000 | for discussions concerning the octupole terms 
of order (ri/r 2 )^). 


It is this perturbative Hamiltonian H' that leads to a secular evolution of the 
triple system. In the hrst approximation, the secular evolution can be described by H' 
doubly-averaged over the mean anomaly of both the inner and the outer orbits. It turns out 
that the doubly-averaged H\ the total angular momentum of the system (ao), the energy 
of each binary (or equivalently oi and 02 ), and the magnitude of the angular momentum of 
the outer binary (/3o, or equivalently 62 ) are conserved quantities. The secular evolution 
of the triple system can therefore be fully described in terms of time-evolving e and gi. A 
useful constant of integration can be derived from equation (H), 


Cll€ 




(Clfi 


- Pir 

4/9o^ 


( 4 ) 


Note that the quantity on the left hand side is conserved as long as the total angular 
momentum (ao) and that of the outer binary (/do) are conserved. Also, whenever i/e -C 2/d, 
Oiecos^J is a constant. 


The secular evolution is a result of coherent additions of perturbations from the tertiary 
to the inner binary. Therefore, any mechanism that perturbs the phase relation of the 
system could modify the evolution of e and gi. For black hole systems, an important effect 
comes from the general relativistic periastron precession of the inner binary. Inclusion of 
this effect can be found in Miller & Hamilton (2002) and Blaes, Lee & Socrates (2002). 

In summary, contribution of the periastron precession in the first-order post-Newtonian 
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approximation can be added to the doubly-averaged Hamiltonian as 
where, 


6 'pn — 


8 X 10“® 


tus/Mq 


^ 2 ^ 1 

aij oi/AU’ 


—A: 6 *pn/ ^/e, 


(5) 


and k = ?)GmQmim 2 a\/(SMib^) is a quantity related to the evolution time scale. Apparently 
the influence of the periastron precession (abbreviated as PN effect) is the largest (in the 
sense that 0 pn is the largest) for systems consisting of a tight and massive inner binary and 
a light, far-out third component. With the addition of this PN effect, the total Hamiltonian 
remains conserved. It can be written as H = kW ([Miller fc Hamilton 2UU2|), where 


^ 9i) = ~2e -f e cos^ / - 1 - 5(1 — e) sin^ 5 fi(cos^ / — 1) -|- 6 'pN/e^'^^, ( 6 ) 


is a conserved quantity. 


2.2. Eccentricity Evolution in Absence of the GR Effect 


The trajectory of ei and gi in the phase space are known to be determined by the values 
of a and (3 for given initial system parameters (jLidov fc Ziglin 1976|) . In the general case 
of initial ei 7 ^ 0,1 and in the absence of the GR effect, the quantities ei, gi undergo cyclic 
oscillations. A necessary condition for a large change in ei is that Iq > cos“^ ~ 39'’ 
for a restricted case where mo -C m 2 -C mi. This condition still holds for arbitrary masses 
when the PN effect is weak (as can be proved using eq. |P). The time scale for the system 
to swing from ei ~ 0 to Ci ~ 1 is estimated to be (jlnnanen, Zheng, Mikkola, fc Valtonen 

rwFj) , 

(MiG (^G _ ,2)3/2 


Tevol ~ 0 . 16 / 


Vm2/ \aiJ (m2/M0)P2 


yr, 


(7) 


where / ~ 0.42 \og{l/ ciq)/\J (sin^(Jo) — 0.4 is a quantity of magnitude a few, Ciq is the 
initial value of ei, and Iq is the initial value of I. This time scale should be much longer 
than the orbital periods of both binaries because it is an orbital averaged effect. 


The minimum value of e (denoted Cmin) for arbitrary masses mg, mi, m 2 , and 
arbitrary eo can be solved for using the fact that the Hamiltonian is conserved, that is, 
hP(eo, 5 'o) = hP(emin,S'™) (eq. 0). We are interested the cases emin ~ 0 and emin'C eo. It 
follows that = 7 r /2 for Jo 7 ^ Jc ( jLidov fc Ziglin 1976|) , e G„ -C eG^, and emin -C 2(3. The 
assumption that emin "C 2(3 is valid as long as there exists a solution for the critical initial 
value cos Jc = —yTo/(2/?). Under these conditions, the cos^ J term in lU(emin,fi'm) can be 
rewritten in terms of cos^ Jo and emin- An explicit approximation formula for emin is then 
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obtained by solving the resulting quadratic equation as 


d/2 


2VL 


6*pn + \/^pN + 20f2eo [cos/q + y^/(2/3)j' 


By using equation (0), we have 


1/2 

^min 


2Q 


+ 


\ 


^PN + 5 


ai (5^ 


( 8 ) 

(9) 


Here 
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n = 5 - 2eo + eo cos^ Jq + + 4eo(cos/o + + 5(1 - eo)(cos^ Jq - 1) sin^(5fo)- (10) 




a/^ 


2/3 


For the restricted case of tuq -C m 2 -C mi, we have /3 ^ 00 , Ic ^ 90°. We recover 
equation (8) of Miller & Hamilton (2002) using their assumptions that Cq ~ 1, 6*pn -C 3 and 
5 cos^ Jo 3. We then also recover the classical limit Cmax = \/3 — C min = ^1 — (5/3) cos^ Jq 
( e.g., Innanen et ah 1997) by setting ^pn = 0. 

It follows from equation (^) that a condition for achieving extremely small Cmin is 
^pN 1 and Jq ~ /c > 90°, which corresponds to a ~ (eq.]^). Note that Jc > 90° means 

that the outer binary should initially be in a retrograde orbit with the inner one. This is 
consistent with what we obtained in the classical limit (0 pn = 0), where a system can evolve 
to e = 1 starting from any initial eccentricity ii a = (3 (or Jq = Jc) (see details in |Lidov fc| 
Ziglin lOTOl) . The critical values of Jq above which the system can have a signihcant change 
in e can also be derived. If we demand Cmin < 0.5eo, we would obtain from equation (P) 
cos^ Jo < 3/5, or Jq > 39° for cq 1 and 6 *pn ~ 0. This is the same as the known classical 
limit for the restricted problem. 


Two important constants can be derived from equation (|^) for triple systems with 
arbitrary masses but negligible PN effect. As ^pn ^ 0, is a constant ~ 3 for Cq ~ 1 and 

h ~ h-, 

5 («2 - 


OiCr 


12 PI • 

It thus follows from equation (^) that the mutual inclination angle at emin (denoted J, 
satishes 


( 11 ) 


cos^ /„ 


3/5, 


( 12 ) 


for ye 

min <2/3. 



















3. Influence of Gravitational Radiation Reaction (GR Effect) 


In this section, we discuss how the GR effect damps the Kozai oscillations and 
eventually leads to the merger of the system. As the gravitational radiation frequency 
depends directly on the quantity oie ( 7 ^, or the semi-latus rectum) (section eq. |^), we 
focus our discussions on the evolution of this quantity. We hrst discuss the equations that 
govern the evolution of the system where the GR effect is included. We then derive the 
minimal e (cmin) obtained at the hrst Kozai cycle and then determine the evolution of Cmin 
during subsequent cycles, depending on the roles of the PN effect and the GR effect. The 
evolution of the frequency of the gravitational wave emitted by the inner binary can then 
be determined. 


3.1. Equations of Motion 

The presence of gravitational radiation has important impacts on the eccentricity 
evolution of the inner binary. First, the radiation carries away energy and angular 
momentum and tends to shrink and circularize the orbit. The decay of the orbit 
subsequently changes the values of a and /3 and thus modifies the secular evolution 
trajectory of the system. In particular, the effect of periastron precession becomes stronger 
rapidly with the decay of the orbit (see eq. [^). This results in rapid oscillations of g\ 
that could destroy the secular oscillation of e. As for the outer binary, the orbit remains 
much wider than the inner one, so its energy (equivalently 02 ) and magnitude of angular 
momentum (/do, or equivalently 62 ) can therefore still be treated as conserved quantities. 

The GR effect is the strongest at e ~ Cmin within each Kozai cycle because of the strong 
dependence of the merger time scale tgr on e (eq. J^). When the PN effect is negligible and 
the GR effect does not affect the Kozai cycles dramatically, we have oi oc according to 
equation (0). The evolution of the parameters QpN-, Tevoi, and tqr with the decay of the 
orbit can then be summarized as follows. 


C\ —4 4 

fc^PN OC oc e,^in 

(13) 

-3/2 3/2 

^evol OC Qi^ OC Cmin 

(14) 

4 7/2 -0.5 

Tgr oc a^e^i^ oc e^;,, . 

(15) 


It is apparent that, as the orbit of the inner binary shrinks and circularizes, the value of tqp 
decreases proportional to while the time the system spends at emin (Rvoi = nevoiA/ejjjijj^) 
increases proportional to ~ e^;^. When tgr becomes of order TevoiA/emin, the gravitational 
merger will take place within one Kozai cycle near e ~ emin- 
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The evolution of the system can be calculated through a set of hybrid equations that 
combines the equations of motion derived from the conserved Hamiltonian W{e,gi) (Lidov 
1976, eq. [32]-[33]) including the PN effect and that from quadrupole gravitational radiation 
( [Peters 1964| , eq. [5.6]-[5.7]). These are the same equations used in Blaes, Lee, & Socrates 
(2002). The evolution equations govern the four parameters e, ( 71 , oi, and ao and are given 
by 


de 

dt 

dgi 

dt 


dai 

dt 

dciQ 

dt 


-lt)KEa\^y/e(l - e)(l - cos^J) sin( 25 (i) + ^^^^(^ - 

6 


304 304 


KEa\^ < —j= 


'a L 


4 cos^ / + (5 cos 2g\ — 1) (e — cos^ /) 

6 *pn 1 


r*os 1 

H—[2 + (1 - e)(3 - 5 COS 2gi)\ + 


/3 


37 


19a?e3-5^6 16^^96^'^’ 

3 Kg 1 /15 7 , Jatii + /3o cos^ / 

- ‘ - -e)- -, 


19 af-5 e2 ^ 8 8 


ao 


where 6 *pn is dehned in equation (^, and 


ke = 7.4554 X (WMe)»" 


Kg = 7.8218 X 10 


Ml' (a2/AU)3(l 

26momiMi^^4 




The evolution of the angular momentum 7 can be written as (see eq. [|T^,[18|), 

d'~1 r 2 • 2 r • o ^ —3 3 —2/ 15 ^ 25 

— = -5KEai sm I sm 2gi - e ~ 8^ ' 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 

( 22 ) 


Equations (|l 6 |) - (p!^ provide a valid description of the evolution as long as the energy 
(oi) is approximately conserved within each cycle, and as the total angular momentum 
vector of the triple system is conserved. These are necessary conditions for the validity of 
the equations derived from the conserved Hamiltonian and the following considerations 
justify them for the systems that interest us. 


We have assumed that the triple systems in globular clusters consist of stellar mass 
BHs with initial orbital separations on the order of AU, and moderate initial e ~ 1 (see the 
discussion of parameters in section | below). The systems are required to obtain extremely 
small emin~ 0 during their secular evolution in order to merge well before the system is 
disrupted by interactions with held stars ([Miller fc Hamilton 2002|) . Equation (^) and (p!7[) 
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indicate that at each Kozai cycle, such systems would spend most of their time at moderate 
e ~ 1 , and a very small fraction of the time (oc ^Cmin) at extremely small values (see also 
Innanen, et al. 1997). Note that gravitational radiation reaction is negligible at moderate e 
and ai as Kg ^ (eq. |^,[^]) for such systems. The energy (equivalently oi) can thus 
be treated approximately as conserved for most of the cycles (see also eq. |T^). 

The total angular momentum vector can also be treated as conserved within each 
cycle as long as emin~ 0. It is known that, for an isolated system under gravitational 
radiation reaction, the angular momentum loss is negligible (oc KG^mfn) e mill 0 (eq. 

second term). Moreover, gravitational radiation reaction will not change the direction of 
the angular momentum of the inner system. It is thus conceivable that as long as emin~ 0, 
the evolution of the angular momentum of the inner binary ( 7 ) is negligibly affected by 
the GR effect, and the vector of the total angular momentum can be treated as conserved 
throughout the cycles. For our systems, the GR effect starts to dominate at emin~ 0 because 
of the extremely small emm required for the system to merge before the system is disrupted 
by a held star. These justify the validity of equations (pTD-(p!^). 



3.2. Estimation of emm 


The minimum value of e that a system can reach within its hrst Kozai cycle (denoted 
emin) is best estimated using equation ([TBI). At e = emin, de/dt = 0 , we obtain. 


sin = 


1 + 


(23) 


where 


7o = 


425 Kg 
3040 ke 


= 1,467 X (“2/AU)’(1 - 4 ?'^ 


niQmiMi \m2) {m2/ 


-AU^’-h 


(24) 


and where sin^ Im is related to emin through equation (|^). We estimate emin by assuming 
that Oi is a constant and solving implicitly the equation 


AW = W{go, eo) - W{gm, Cmin) = 0. (25) 

In our numerical work, we actually solve for emin by minimizing (AhF)^ using a downhill 
simplex method developed by Nelder and Mead (1965) (available in Matlab). This method 
has proved to be robust in hnding local minima. However, because of the existence of 
multiple solutions to equation ( p5[) (with at least one other solution e = eo and gi = go), a 
good initial guess is essential for a well-converged solution. 









11 


We made our initial guess emin by considering the following two limits for based on 
equation (^). In the hrst limit, gravitational radiation reaction is weak, or S> 70 , 

we have gm ~ 90°. Equation (H) can then serve as a very good approximation to emin- In 
the second limit, gravitational radiation is very important, gm could deviate signihcantly 
from 90°. The extreme limit occurs around Jq = Ic or = Pq (eq. 0). As with the critical 
case of a = /9 in the classical limit, a system starting with cq ~ 1 will pass through an 
unstable stationary point e ~ 1 and gi ~ sin“^ ^2/5 and move on towards a stationary 

(and stable in the classical limit) point e ~ 0 , and gc = sin“^ 


1/(2/5)y'( 8/32 - l)/5j ( iLidov] 
|fc Ziglin 1976| , eq. [46]). If it were not for the GR effect, the system would reach ei = 1 
at this stable point. With the presence of the GR effect, the angle gi will be steered away 
from this point while e is steered away from a true zero value (see also Fig. However, 
we expect the value of gm to be very close to gc (at emin~ 0 , dg/dt oc e-mm ~ 0 ) (eq. |1T7| , 
ignoring the PN effect). In this limit, the minimum e can thus be estimated using gm ~ gc 
and assuming constant oi. Specihcally, we again use the minimization program introduced 
previously to solve equation 


^min bo ■ 2 T • o ’ 

sm ImSm2gc 

for Cinin, where sin^ Im is related to by equation (H). 

Without knowing the solution ahead of time, we simply feed values of emin estimated for 
these two limits as initial guesses to the minimization program and then select the solutions 
of {AWy = 0 with minimum residues. In all cases, the residues of the minimization are 
monitored to ensure correct convergence. 



3.3. PN effect 


The importance of the PN effect increases rapidly with the decay of the orbit as 
6 *pn oc ah^. This effect has a direct impact on the oscillatory behavior of gi (eq. 0 ), 
through which it can affect the evolution of e and thus aiCmin (eq. [HI]). The importance 
of the PN effect to the evolution of e and aiCmirn however, depends on the role of the GR 
effect. 


The PN effect can be neglected if the GR effect dominates in the evolution of e before 
the PN effect becomes important in the evolution of gi. Once the GR effect dominates, the 
evolution of e decouples from that of gi and the Kozai cycle is terminated. In the cases 
we are interested in, this transition occurs near emin~ 0. When both the GR and the PN 
effects are not important, the quantity aiCmin is roughly a constant (eq. [|TT|). After the GR 
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effect dominates, it evolves according to equation (5.11) of Peters (1964) and remains a 
constant as long as einin~ 0. The overall evolution of Oiemin can therefore be conveniently 
described using Peters’ equation with oi and emin of the first Kozai cycle (see section ^.21 ) 
as the initial data. 


The PN effect has a signihcant impact on the evolution of e and aiCmin if it dominates 
in the evolution of gi before the GR effect becomes important. In this case, a dominating 
PN effect will introduce a fast oscillation in gi and thus in e (eq. 0) and will terminate 
the Kozai cycle. The GR effect will eventually dominate in the evolution of e and Peters’ 
formula can again describe the evolution of Oiemin- However, when the PN effect is 
important, the quantity aiemin can no longer be approximated as a constant before the GR 
effect dominates. It will evolve to a larger value than that in the hrst Kozai cycle (see 
eq. ||3.2|| ). The results based on Peters’ formula with initial values taken at the hrst Kozai 
cycle can therefore be used to set the lower bound on aiemin- 


An upper limit for the quantity aie^i^ when the PN effect is important can be set based 
on the following considerations. First, the fast oscillation induced by the PN effect will be 
damped away under the GR effect. In light of equation (^, the quantity G will evolve from 
a rough constant of 3-5 towards G ~ 6*pn/a/^ — cx) as 6*pn —>■ oo with the decay of the 
orbit. As a result, emin —^o, that is, the oscillation will be quenched. The time scale for the 
system under the PN effect to swing from ^fi ~ 0 to 7r/2 with a small change in e and Oi is 


TpN ~ 6 /(kp;0pnU^ ^) ~ 5 X 10^ 


VAUy (Ml/Mo) 1-5 


yr. 


(27) 


Secondly, neglecting the GR effect, the evolutions of /, and gi under a dominating PN effect 
can be written as 


— = —5Kp;a^sin/cos/(aie) ^'^^sin 25 fi, 

dgi I 5 6 *pn 

— = Keap -. 

at e 


(28) 

(29) 


where we have applied equation (^ and (^^ to obtain equation (|2^) . It is apparent that 
the oscillation of / with gi is also damped rapidly as 6 'pn/v^ —oo with the decay of the 
orbit. Integrating over gi = [0,7r/2] for dl/dgi, we obtained an estimation for the oscillation 
amplitude 


SI ^ 


2.5^sin2/. 


7PN 


(30) 


An upper limit of SI can be estimated by the fact that ^ > 10 cos^ ~ 6 as required 
for a dominating PN effect (eq. 0.0). and that sin 2/ < 1. Here, we have used the fact 
that the PN effect becomes dominant near emin (or J^) where oi decreases most signihcantly 
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under the GR effect. We thus estimate that under the dominating PN effect, we have 
Im < 39° + max(5J) < 63°. An upper limit for the quantity can therefore be set 

according to equation (^) 

5 




PI 


(31) 


By comparing this limit with equation (|^, we simply demand that the upper limit of 
Oiemin be three times its lower limit obtained with Peters’ formula. 


For each given set of initial system parameters, we consider the PN effect to be 
important if aiCmin obtained within the first Kozai cycle, which set the lower limit of this 
quantity during its evolution, satisfies 

> 7o , (32) 

where 6*pj,j = is a constant. This condition is obtained by the requirement that the 

GR effect is not important for the evolution of e (eq. 0) when the PN effect starts to 
dominate the evolution of g\ (eq. 0), that is, 

> 7o, (33) 

^PN ^ 6\/e. (34) 


Equation ( |3^ and the comparison of tpn (eq. [0) with tqr (eq. |^) indicates that the PN 
effect is important only for systems with sufficiently large Oi and emin- As a result, the PN 
effect is relevant only in low eccentricity cases in the LIGO band (cf. eq. |l37|| ). 


3.4. Evolution of Gravitational Wave Frequency 


The frequency of gravitational waves (/gw) emitted from an inspiraling binary depends 
strongly on its orbital frequency /orb and eccentricity Ci, where 

_ VGMi 

/orb- • 

In the limit of the quadrupole approximation, the power at the nth harmonic of /orb have 
been derived in equations (19) and (20) of Peters & Mathews (1963) in terms of spherical 
Bessel functions. For a circular orbit, the power is concentrated in the second harmonic. 
For an eccentric orbit, the power spreads over a broad frequency band and the maximum 
occurs at a harmonic rim S> 1 for e ~ 0. For instance, this peak frequency could occur at 
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the n-m = 10® harmonic of the orbital frequency for e = 2 x 10 In this section, we discuss 
the evolution of the peak frequency with the eccentricity e\ (or e). 

We have derived an approximate expression for the “peak harmonic” of order rim in 
terms of e for convenience of our calculation as follows. We hrst used the formulae in Peters 
& Mathews (1963) to calculate Um for various values of e = 10“®-1 and we then applied a 
least square £t to the data. An excellent £t to ^^(e) was obtained, 


— 


2(1 + Cl) 


1.1954 


(l-ef)i- 


(36) 


This formula preserves the correct limit of = 2 for the circular case of Ci = 0. The 


relation of the peak frequency to e is obtained by combining equations (|35|) and (^6]), 




TT 


(o,e)i^5 


(37) 


It indicates that the peak frequency i^ inversely proportional to the period of a circular 
Keplerian orbit with radius equal to the actual orbit’s semi-latus rectum oie. 

We actually are interested in ci as a function of rather than /^(ei), since the 
focus of this article is the orbital eccentricity as the GW frequency evolves through the 
LIGO band. We compute ei from equation (0. by applying our knowledge of the evolution 
of aie. Let the initial system parameters be ai = oiq, e = cq (ei = ciq), and the minimum 
e expected for the hrst Kozai cycle be (eq. |^) and Cm = y^l — e ® J . We consider 
the following three possibilities separately, (i) The frequency enters the LIGO frequency 
band before e reaches In this case, oi ~ oio, and e lies in the range (ii) The 

frequency enters the LIGO band after was reached and the GR effects dominates in the 
evolution of e before the PN effect becomes relevant. The evolution of aie is then the same 
as for an isolated system evolving under gravitational radiation reaction with initial values 
of Oi = Oio and e = emm- That is, (eq. [5.11] in Peters 1964), 


0 W ei y2/i9 / 1 + i21/304e? 
(«ioemin) j (l + 121/304e^ 


870/2299 


(38) 


(iii) The frequency enters the LIGO band after was reached, but the PN effect becomes 
important before the GR effect dominates in the evolution of e. In this case, the lower limit 
of aie is obtained from equation (^) and the upper limit is set to be three times this lower 
limit. The bounds for ei at a hxed /gw can be obtained by applying these two limits. 
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3.5. Numerical Examples 


The numerical evolutions of e and I are shown for two typical examples in Fig. ^ 
and in Fig. ^ for illustration. These evolutions were obtained by integrating the ordinary 
differentiation equations ('[l 6 [)-(|l^) for given initial data and then using equation (||) to 
derive the mutual inclination angle I. We used a medium order Runge-Kutta method 
(available in Matlab as ode45) for the integration. 


The first example (Fig. represents a case in which the system merges after many 
Kozai cycles. This is also a case where the PN effect becomes important before the 
GR effect dominates in the evolution of e. The system parameters are chosen to be 
mo = mi = 1712 = IOMq, Ci = 62 = 0.1, Oi = 10 AU, 02/01 = 10, gi = 0, and Iq = 95.3°. 
The gradual increase in the time scale of the Kozai cycle is apparent. For most of the 
cycles, Cmax ~ 1 (emin ~ 0). Near the emin of the last Kozai cycle, fast damped oscillations 
of e (and I) due to the PN periastron precession are visible before a monotonically increase 
of e due to the dominating GR effect (see also Fig. H). The value of / experiences damped 
oscillation after the PN effect becomes important and converges to a constant after the GR 
effect dominates in the evolution. Note also that before the last cycle, the minimum I is 
around 39° as predicted by equation ([T^). 


We show in Fig. for the hrst example, the evolution of e with oi (solid line) and 
compare with the results based on the lower limit of aiemin obtained with equation 


(dashed 


and a upper limit that is three times this lower limit as discussed in section 
lines). Also shown in dashed-dot lines are the evolutions of the time scales of the GR and 
the PN effects normalized by the Kozai time scale The evolution path of e starts from 
the upper right of the solid line and moves towards the left. It is apparent that oi remains 
essentially a constant within each cycle before TcR/Tg^oi ~ 1 - The decay of oi occurs mainly 
near Cmin of each cycle. Damped oscillations are apparent near rpN/r/^oi ~ 1- There is an 
excellent agreement between our numerical results and the prediction in sections 


The second example (Fig. ^ represents a case where the system merges within the hrst 
Kozai cycle. The system parameters are chosen to be mo = mi = m 2 = 5Mq, ci = 0.01, 
62 = 0.51, Oi = 2.1920 AU, 02/01 = 5, gi = 0, and Iq = 99.2°. This is a case similar to 
the classical a = case. That is, the system evolves along e ~ 1 until it passes through 
an unstable stationary point at gi ~ sin“^ \/2/5 ~ 39° and then evolves towards extremely 
small values of e (Fig. solid line with open circles). The evolution is then quickly 
dominated by the gravitational radiation reaction near Cmin- 


The evolution of e with oi (solid circles) for this second example is shown in Fig. p. 
The predicted e min of the hrst Kozai cycle calculated with the method in section at 
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initial a\ = Qiq is shown in symbol ‘x’. The expected evolution of e with ai by equation (^81) 
starting from this predicted emin and Oiq is shown by the dashed line. It is clear that there is 
an excellent agreement between the predicted evolution of aiCmin and the numerical values. 
The same excellent agreement has been found for all such cases we have investigated. 

The phase diagram for the evolution path of e vs gi for systems with the parameters of 
our second example but with different initial values Jq are shown in Fig. It is apparent 
that as long as Jq is away from the critical values Jc, the e^am is reached around gi ~ 90". 
Near the critical value, gm ~ gc- 

The evolutions of eccentricities ei as a function of gravitational wave frequencies 
for our hrst and second examples are shown in Fig. and (||). In both cases, the GW 
frequency spans a wide range of eight orders of magnitude. This means that the GW’s from 
these sources can be good candidates for detections in the LISA band as well as LIGO. 
For the case of many Kozai cycles in Fig. |^, the eccentricity drops sharply from ei ~ 1 at 
/gw ~ 0.1 Hz (the end of the LISA band) where the GR effect starts to dominate, and is 
nearly zero at 10 Hz (the beginning of the LIGO band). This is expected as equations (|T7D 
and (pH]) predict that the eccentricities drop with the frequency roughly proportional to 
/'"gw extreme cases as shown in Fig. j^, the eccentricity remains at the signihcant 

value of ~ 0.9 at /q^ = 10 Hz and could be extremely high at lower frequencies. 


4. Eccentricity Distribution 


The Kozai mechanism can drive the inner binaries of triple systems to extremely small 
e and merge before disruption by interactions with held stars. The time scale for disruption 
(the same as the stellar encounter time scale) is given as (jMiller fc Hamilton 2002]) 


„ . p,5 _iAU IOMq 

Tenc ~ 6 X lO^’u/-yr. 


a.2 M 2 


(39) 


where the number of stars in the globular cluster \s N = lO^ne. We assume ng = 1 in this 
paper. Successful mergers of this sort therefore require 


"^evol ^ ^enc) 


'rGR(Ol; 


Y^mir 


< Tpi 


(40) 

(41) 


Equation (|40|) ensures that the system has completed at least a half Kozai cycle to reach 
small e. Equation (Plj ) ensures that the system can merge successfully before the next 
interruption. The term on the left-hand side of equation ( 0 ) represents the total lifetime of 
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the system which is roughly a factor of the merger time scale tgr. This is because 

the system only spends a signihcant fraction of yTTTn hs time near emin where the GR effect 
is the strongest. 


4.1. Distribution of Initial e min and ai 


We have explored the parameter space of emin and oi for successful mergers driven by 
the Kozai mechanism in globular clusters restricting ourselves to a triple system consisting 
of three 10 Mq BHs. We £x the initial parameters gi = 0, and Cig = 620 = 0.01, and we 
have chosen 02/01 = [20,10,5,3] and assume a uniform distribution in the initial mutual 
inclination angles Jq and in the semimajor axis oi. The minimum initial value of oi was set 
to 0.2 AU so that during its interaction with a third body, the recoil velocity associated 
with binary hardening is low enough for the system to remain in the cluster ([Miller fc 
[Hamilton 2002|) . There was no presumed upper limit for Oi. For each given oi, Jq, and 
02 / 01 , we hrst calculated emin with equation (^5[) for each given oi and Jq. Equations (|40|) 
and dl) were then evaluated to determine the permitted parameter space. 


We show in Fig. || the permitted values for ai (upper set) and corresponding Cmm 
(lower set) vs Jq. There are four sets of cone-shaped distributions, each corresponding to a 
given 02 / 01 , with value decreasing from the left to the right. Each distribution centers on 
the critical value Iq ~ Ic which depends solely on 02/01 (eq. @]) as cq is hxed. There is a 
cut-off at a maximum Oi. This cut-off in oi is set by the constraint from equation ([40[) . For 
a hxed O 2 / 01 , larger oi implies larger 02 and therefore shorter Tenc (eq. f^)- Beyond the 
maximum Oi, the time scale Tenc becomes so short that the system is disrupted before it can 
hnish a half Kozai cycle to reach emin- For the same reason, a larger 02/01 leads to smaller 
cut-off in oi. 


The overall cone shape in oi and emin results from the fact that the merger time scale 
should be less than the encounter time scale (eq. |^^). Near the center of the cone where 
-fo ~ Ici the system can reach extremely small Cmin (eq. [§]) which allows maximum possible 
values of oi under the constraint of equation (pj). The further Iq is away from the center 
{Ic), the larger Cmin is, and so the smaller are the Oi values that satisfy equation (^). 


The range of Cmin within each distribution is caused by the PN effect. In the classical 
limit (with no GR or PN effect), emin depends only on the values of Iq and a^jax. It is 
independent of oi for each distribution, making emin have a unique value at each Iq. This 
is apparent for e at Iq far away from center of each distribution for 02/01 = 3,5,10 in 
Fig. ||. A significant contribution from the PN effect makes emin depend on oi, giving emin 
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a finite range (eq. p,[Q). The PN effect is especially pronounced near the center of each 
cone (Jq ~ Ic), where the contribution from the classical Hamiltonian is relatively smaller 
(eq. [§]). The finite range of Cmin is more pronounced for cones with larger 02/01 values as 
the PN effect is more pronounced at larger 02 / 01 . This also explains the overall large emin 
at larger 02 / 01 . 

The ejnin values estimated including the GR effect can be several orders of magnitude 
larger than those without the GR effect in the region Jq ~ Jc, where the GR effect is the 
strongest. This makes little difference in finding the permitted parameter space for oi, as 
6 min will be extremely small in this region and equations (^) will be satished in either 
case. However, the evolution of gravitational wave frequency depends sensitively on our 
knowledge of Oiemin- R is therefore necessary to include the GR effect. 


4.2. Distribution of Eccentricities in the LIGO band 


We proceed to investigate the distribution of eccentricities when the emitted GW 
wave enters the LIGO band. We restrict attention to stellar mass black holes in triple 
systems and show only a representative case of triple systems consisting of three 10 Mq 
black holes. We consider a uniform distribution of initial mutual inclination angle Jq, initial 
eccentricities of the inner and outer binaries, ratios of the semimajor axis 02 / 01 , and initial 
gi- A summary of the parameter space we have investigated and the representative values 
we have included can be found in Table |I|. 

We choose the upper limit of O 2/01 to be 30 based on equations (|^, as the parameter 
range for oi diminishes with larger 02/01 (see Fig. ^). A general upper limit for the ratio 
02/01 can also be set based on the requirement that eo > Cmin; h follows from equation (|) 
that 



(42) 


This expression is similar to the limit given by Blaes, Lee & Socrates (2002) (eq. [5]) with 
a slight difference in the constant term and difference in the power of e ( 1/6 here vs 1/2 
there). The lower limit for 02/01 is set to be 3 based on the fact that the triple system is 
stable only if ([Mardling fc Aarseth 2001]) , 



(43) 


We follow the same procedure described in the previous subsection. We first calculate 
the permitted parameter space for emin and oi. The results (not shown) are many versions 
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of the same distributions of those shown in Fig. ^for various a^jdx. The span of permissible 
Jo is roughly 90-106°, as required to obtain extremely small emin- We then calculate the 
expected eccentricity e\ or its upper and lower limit at given frequency following the 
numerical procedures described in section taking into account of both the GR and PN 
effects. 


The number distribution of eccentricities at ~ 10 Hz is shown as a histogram 
in the left panel of Fig. |^. The distribution is shown in percentage relative to the total 
number of systems that would merge before interruption by field stars. The histogram in 
grey shade represents the distribution without considering the impact of the PN effect on 
the evolution of aie-^am (see discussions in section ^.3| ). The histogram in dark shade includes 
the upper limit due to the impact of the PN effect. The x-axes of the histograms is shifted 
for a better view. It is apparent that the PN effect affects mostly the low-eccentricity 
systems and it makes little difference in the histogram. 

At /qyv = 10 around 70% of the merger systems have eccentricity < 0.1, a little 
more than 50% have eccentricities < 0.05. About ~ 2 % have eccentricities ei ~ 1, most of 
which are those with initial Jq very close to the critical angle Jc, as the GW’s reach 10 Hz 
before the systems reach emin within the first Kozai cycle. These are the type of systems 
studied in Fig. I, |,i and ||. 

We also show eccentricity distributions in the right panel of Fig. |0|for/-w = 40 Hz 
(grey shade) and /q^ = 200 Hz (dark shade). At 40 Hz, all eccentricities are well below 0.2. 
At 200 Hz, all are well below 0.02. This is not surprising as equations (^7]) and (|38D predict 
that the eccentricities drop with the frequency roughly proportional to This is 

consistent with the fact that the GR effect circularizes the system in a very rapid fashion. 


5. Discussion 

We have studied the evolution of eccentricities of black hole mergers driven by the Kozai 
mechanism in triple systems residing within globular clusters. Eccentricity distributions 
at gravitational wave frequencies relevant to LIGO have been presented. The evolution 
of these systems were investigated including the Kozai mechanism, the post-Newtonian 
periastron precession effect, and gravitational radiation reaction. 

We conclude that around 30% of the systems possess eccentricities > 0.1 when the 
emitted gravitational waves reach 10 Hz. Around 2% of the systems are extremely eccentric 
at 10 Hz. However almost all our merger systems possess eccentricities well below 0.2 at 
40 Hz, and below 0.02 at 200 Hz (see also a brief discussion in Miller 2002). These merger 
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systems, on the other hand, promise extremely high eccentricities at the lower frequencies 
of the LISA band (see Fig. |0],[il)- 


About 30% of our merger systems possess e > 0.1 at 10 Hz, the lower end of the 
advanced LIGO frequency band ([Fritschel 2002|) . It is thus important to determine whether 
it is possible to use LIGO’s current circular-binary search templates to detect such systems. 
For highly eccentric systems, the higher harmonics probably need to be taken into account 
when optimizing the search templates. It is plausible that the eccentricities associated 
with these systems are not important at 200 Hz or perhaps even at 40 Hz. It is important 
however, to determine the limit of ci below which, circular-binary templates need to be 
replaced by new, eccentric-binary templates (see [Martel fc Poisson 1999|). 


The eccentricity distribution was calculated for triple systems consisting of three 10 
Mq black holes. Similar eccentricity distribution, however, can be found for triple systems 
with individual masses in the range of 3-25 Mq known for the observed galactic stellar 
mass BHs ([Bailyn, Jain, Goppi, fc Grosz 1998]) . We have assumed that, in globular clusters, 
triple systems with one or more BHs of mass > 25Mq are rare. The mass parameters affect 
very little the values of Cmin a system can reach (eq. p],p,|Q). Its effect to the shape 
of eccentricity distribution is also very weak compared with other parameters such as oi, 
02 / 01 , and e (eq. |^, |^). Any effects due to different masses will be further averaged out 
if we include a uniform distribution of masses within this range. 


Our calculation of the gravitational radiation reaction is based on the Newtonian 
quadrupole approximation. At 10 Hz, this approximation is still valid as n/c ~ 0.1 at the 
periastron. However, at higher frequencies, especially near 200 Hz in the LIGO frequency 
band, this approximation breaks down as the speed of the binary orbit is approaching that 
of light. However, all our binaries become so circular well before they reach 200 Hz, that 
our results are probably still relevant. 
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Fig. 1.— Geometry of a hierarchical triple system. 








24 




0 5 10 15 

Time (x 1 yr) 


Fig. 2.— Secular evolution of e = 1 — of the inner binary and mutual inclination 
angle I between the inner and outer binaries. This evolution was computed by integrating 
equations ([I6|) -([l9| ) which include contributions from Kozai mechanism, PN periastron 
precession, and gravitational radiation reaction. The initial system parameters are: mo = 
mi = m 2 = 10 Mq, Cl = 62 = 0.1, ai = 10 AU, 02/01 = 10, gi = 0, and Jq = 95.3°. This is 
a typical case in which the PN effect becomes important before the GR effect dominates in 
the evolution of e, and that the system merges after many Kozai cycles. 
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(AU) 


Fig. 3.— Evolution of e (solid line) and the time scales of gravitational merger and the 
PN effect (dashed-dotted lines) with the semimajor axis oi for the same initial parameters 
and evolution equations as in Fig. ^ The evolution path of e starts from the upper right of 
the solid line and moves towards the left. The upper and lower limits on emin(ai) based on 
equations m and (|3l| ) (dashed lines) are also shown. 
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Fig. 4.— Secular evolution for e of the inner binary and of /, as computed by integrating the 
evolution equations (0)-® which include contributions from the Kozai mechanism, PN 
periastron precession, and gravitational radiation reaction. The initial system parameters 
are: mo = mi = m2 = 5 Mq, ei = 0.01,62 = 0.51, oi = 2.1920 AU, 02/01 = 5, 5^1 = 0, and 
Jo = 99.2°. This is a typical case where the GR effect dominates in the evolution within one 
Kozai cycle. 
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Fig. 5.— Phase diagram of e vs gi for the same initial parameters as in Fig. ^ but with 
various Iq. We have Iq = 98° (dashed line), Iq = 99.05° (dotted line), Jq = 99.2° (solid line 
with circles), Iq = 99.3° (dash-dot line), Iq = 99.5° (solid line). For general cases, Cmin occurs 
at Qi = Qm = 90°. For Iq ~ Jc, gm is significantly deviated from 90°. 
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Fig. 6.— Evolution of e with the semimajor axis ai (open circles) for the same initial 
parameters and evolution equations as in Fig. |^. Note that the evolution path starts from 
the upper right and moves towards the left. The predicted Cmin by equation (1^) at Oig is 
shown by the symbol ‘X’. The predicted evolution of e vs Oi from equations (l38|) is shown 
by the dashed line. There is an excellent agreement between the numerical values and the 
predictions. This is a typical case where the GR effect dominates within one Kozai cycle 
and the PN effect is negligible. 
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Fig. 7.— Evolution of ei with gravitational wave frequency for the system shown in 
Fig. i The spans eight orders of magnitude. The eccentricity is nearly zero at and 
above 10 Hz. 
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Fig. 8.— Evolution of ei with for the system shown in Fig. The again spans 
eight orders of magnitude. The eccentricity is around 0.9 at 10 Hz. 




Fig. 9.— Parameter space for ai (upper set) and corresponding values of emin (lower set) vs 
Jo for 02/01 = 20,10, 5, 3 from left to right (see text). 
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Fig. 10.— The distribution of eccentricities at = 10 (l^ft panel) and 40 Hz, 200 
Hz (right panel) for black hole mergers driven by the Kozai effect in a globular cluster. See 
Table |l| for the parameter ranges used. Each parameter is uniformly distributed in its range 
unless specihed. 
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Table 1. Parameter ranges investigated for Kozai-mechanism-driven BH mergers in 

globular cluster 


Mbu (Mq) 

ai (AU) 

02/01 

eo gi{°) 

Io{°) na 

10 

0.2-30 

3,5,10,20,30 

0.01-0.901 0-90 

85-110 1 


Note. — Mbh refers to black hole masses for all three components; cq 
refers to eccentricities for both inner and outer binaries. Numbers of steps 
for ai, eo, gi, and Jq are 60, 4, 4, and 100 respectively. See the text for 
definitions. 



